Theory of Coherent and Incoherent Nuclear Spin-Dephasing in the Heart 
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We present an analytical theory of susceptibility induced nuclear spin dephasing in the capillary 
network of myocardium. Using a strong collision approach, equations are obtained for the relaxation 
rate of the free induction and the spin echo decay. Simulation and experimental data are well 
predicted by the theory. Since paramagnetic deoxyhemoglobin as the origin of nuclear spin dephasing 
has a higher tissue concentration in myocardium supplied by a stenotic, i. e. significantly narrowed, 
coronary artery, spin dephasing might serve as a diagnostic tool. Our approach can be modified for 
capillary networks in other tissues than myocardium and may be applied in material science. 

87.59.Pw, 76.60.Jx 



In cardiology a paramount goal is the detection of sig- 
nificant narrowing, called stenosis, of coronary arteries, 
and - even more important - the detection of the asso- 
ciated area of the cardiac muscle. There is evidence that 
this area can be distinguished from normally supplied 
tissue without contrast agent by NMR imaging from its 
different transverse spin relaxation time H]. Water pro- 
ton spin dephasing is induced by spin-spin interactions 
and susceptibility differences between different histologi- 
cal substructures such as intra- and extravascular space. 
In native tissue this difference is related to the para- 
magnetic property of deoxygenated hemoglobin which 
induces perivascular field gradients. Deoxyhemoglobin 
which results from intracapillary dexoygenation of al- 
most 100% oxygenated arterial hemoglobin, is predom- 
inately located in the capillary and venous system. In 
the cardiac muscle this blood oxygen level dependent 
(BOLD) effect is almost solely due to the magnetic field 
inhomogenties around capillaries since they contain more 
than 90% of intramyocardial blood volume |Q. In my- 
ocardium supplied by a stenotic coronary artery, coro- 
nary autoregulation tries to maintain sufficient perfusion 
by conrpensatory filling of capillaries (recruitment), i.e. 
deoxyhemoglobin increases in tissue and spin dephasing 
is accelerated when compared to normally supplied tissue 
PI . This difference of spin dephasing can be pronounced 
by pharmacological coronary dilation since a higher oxy- 
gen supply decreases the tissue concentration of deoxy- 
hemoglobin in normal myocardium P-p| whereas it does 
not in the poststenotic myocardium observed ||l|. 

Though the relation between spin dephasing and coro- 
nary microcirculation is of great interest for the diagnosis 
of coronary artery disease, at present there is no satisfy- 
ing analytical theory that describes the dephasing of dif- 
fusing spins in the inhomoegneous magnetic field of the 
capillary system. There exist theories on spin dephasing 



for very long and very short correlation times of diffusion 
induced magnetic field fluctuations. The assumption of 
static dephasing of spins as proposed by Yablonski et al. 
B is not valid in human cardiac NMR imaging systems 
(field strength less than 2 Tesla) since the diffusion dis- 
tance I during relaxation (relaxation time Tj* > 30ms, 
D — 1/im-^/ms, I > yj2T^D k. 1 [iva) has the same mag- 
nitude as the intercapillary distance (19/im). Thus, a 
single nuclear spin experiences almost the whole pericap- 
illary magnetic field gradient during its relaxation. Kisse- 
lev and Posse [0 recently proposed analytical treatments 
for long and short correlation tinres. The first extends 
the static dephasing model by including diffusion. How- 
ever, only the effect of a linear local field gradient on 
relaxation is considered. Application of this approach 
would be valid if large vessels would dominate the mi- 
crovascular architecture of myocardium, which is, how- 
ever, not the case. The other is a perturbation approach 
in the local fields. In the cardiac muscle the motional nar- 
rowing condition r(cj^)^" <C 1 with correlation time r 
and variance of fiuctuations of the local fields (expressed 
as frequency) [uP') is not fulfilled, though. Approaches 
which consider spin dephasing for arbitrary correlation 
times as the Gaussian or Lorentzian phase approxinra- 
tion were not sucessful since they neither revealed the 
correct field dependence of the relaxation rate |q], nor 
the correct fiuctuation-dissipation relation W . 

In this paper an analytical theory is proposed which 
describes coherent and incoherent spin dephasing in the 
capillary network of myocardium. It is not an extension 
of the static dephasing model or a perturbation approach. 
Instead, it is based on a strong collision approximation 
which includes both, the static dephasing and the mo- 
tional narrowing regime as two limiting cases. 

It has to be emphasized that this approach is not re- 
stricted to myocardium, but can be modified for the de- 



scription of transverse relaxation in capillary networks of 
other tissues as well. Applications in material science are 
possible, too. For example, the effect of complex poten- 
tials on spin diffusion, and, hence, on spin dephasing, is 
exploited to determine the niicrostructure of spatially or 
magnetically heterogeneous media H . Examples for the 
former are porous structures, for the latter polymer solu- 
tions and liquid crystalls in which magnetic heterogenity 
is induced by either paramagnetic inpurities or an in- 
homogeneous electrical current. After modifications for 
the appropriate geometry the analysis presented below 
should also be applicable in that field of research. 

We will first discuss the tissue model employed. The 
regular parallel architecture of muscle fibers and capil- 
laries in myocardium ||l^ allows to reduce considerations 
from the whole tissue to one tubular capillary (radius 
Re) which is concentrically surrounded by its mean cylin- 
drical supply region with radius Rs (Krogh model Q). 
Nuclear spin diffusion in the whole tissue is replaced by 
restricted diffusion in the supply region, i.e. reflectory 
boundary conditions are introduced at Rs, and, since 
the capillary wall is nearly impermeable for water on the 
time scale of spin dephasing ||l2[, also at Re- We only 
consider relaxation of extracapillary nuclear spins since 
in myocardium the relative intracapillary blood fraction 
<^ = Rl/Rl is less than 10%. The effect of diffusion in 
field gradients along the capillary axis on spin dephasing 
is negligible, i.e. it is only necessary to consider two- 
dimensional diffusion of spins within two concentric cir- 
cles {Re, Rs)- Basic magnetostatics provides the suscep- 
tibility related spin frequency shift around the capillary 
in cylindrical coordinates as 



uj{r,<l)) = -~5ujRlcos{2(j))/r^ 



(1) 



with the characteristic equatorial frequency shift 5lj = 
2-11^ AxBq cos{9) where 7 is the gyromagnetic ratio. Ax 
the susceptibilty difference and 9 is the tilt angle between 
the capillary axis and the external field which is almost 
90° in a clinical scanner. SinceAx < 8 ■ 10~* for blood 
the intracapillary magnetization is less than 1 mgauss for 
clinical NMR magnets and 5u! < 168 rads~^. The con- 
tribution of fields of neighbouring capillaries is negligible 
due to the intercapillary distance of 19/im. 

The starting point of our analytical treatment is the 
free induction decay of the local nuclear transverse mag- 
netization. The time evolution (polar notation: m = 
TUx — iruy) ai X — (r, 0) is determined according to the 
Bloch Torrey diffusion equations p3[ which, written in 
rotating frame coordinates, are 



dtm{x, t) = {DV^ + iLu{x)) m{x, t) 



(2) 



with diffusion coefficient D. The formal solution of 
Eq. (||) is m{x,t) = exp[(DV2 + iuj{x))t]m{x,Q). Since 
the spatial resolution of clinical NMR scanners is far 
above that of capillary dimensions (mm vs ^m) the 



observable is the mean transverse magnetization M{t) 
within the supply region with volume V 



M{t) = V^^ I dx exp[(£iV2 
Jv 



iuj{x))t\m{x,{)) . (3) 



The homogeneous water concentration in the cardiac 
muscle implies a homogeneous initial magnetization dis- 
tribution m{x, 0) which for simplicity is normalized to 
r7i(a;,0) — 1. Since Eq. (||)cannot be solved analytically 
when the local frequency is given by Eq. (Q), a different 
approach is necessary. 

Let us first consider the fluctuations of the local mag- 
netic fields, represented by the corresponding local NMR- 
frequencies lij{x), that a diffusing nuclear spin experi- 
ences. Their autocorrelation function K{t) is deter- 
mined by K{t) — V^^ Jy dx uj{x) exp{tDV^)uj{x) , where 
exp(fZ3V^) is the Green's function of the restricted dif- 
fusion. The correlation time r of the field fluctuations 
is defined as the mean relaxation time fl^ of the auto- 
correlation function, ''" = /„ dtK (t) / K (0) . Using the 
techniques of Rcf . [|lj| , we obtain ||l^ 



= {Rl/AD)H,)/{,-l) 



(4) 



Inserting realistic values of relative intracapillar blood 
volume (<j = 5 to 10%), capillary radius {Re = 2.5 to 
2.75/im) and the diffusion coefficient {D — 1/xm^/ms) 
into Eq. (Q) implies that r < 6ms. Relaxation times of 
the free induction decay observed using clinical scanners 
are Tj* > 30ms. Hence, the relaxation time describing 
the contribution of field inhomogenities is even longer. 
This implies that on a time scale in which significant 
variations of M{t) occur the local field fiuctuations are 
stochastically independent. 

This stochastic independence suggests replacing the 
diffusion operator DV^ in Eq. (H) by a strong collision 
operator mm. This operator describes the diffusion as a 
stationary Markov process with a single rate A that gov- 
erns the relaxation to the steady state distribution p{x), 
independent of the initial state. Due to the homogeneity 
of water in tissue the steady state distribution is sim- 
ply p{x) = l/V. The generator D of the corresponding 
Markov process is then 



D = A(n - id) 



(5) 



where II denotes the projection operator onto the func- 
tional space generated by p{x), i.e. the application 
of this operator on some function f{x) is, 11/ (x) = 
V~^ Jydx f{x) . The constant A is determined selfcon- 
sistently by requiring that D describes correctly the field 
fluctuations. This implies A = t~^ , see Eq. (H). 

Instead of directly solving the time evolution of the 
magnetization, Eq. (W, under the strong collision as- 
sumption, it is more convenient to consider the Laplace 
transform 



M{s)^V-^ I dx [s - DV"^ - iuj{x)] ^ m{x,0) . (6) 
Jv 

Replacing DV^ by D and utilizing the operator identity 
{A + By^ = A^^-A^^B{A + By^ we obtain after some 
algebra 



M(s) = (l + V(s+^^F + ^^^+ 



?V(s + T-i)2 + 5cj2-r-i(l + 



(7) 



M{t) can be obtained from this result to arbritrary 
accuracy as a multi-exponential function M{t) w 
X]iy=i fi^ exp(— r,yt), the parameters determined from (|7|) 
using the generalized moment method of p4[ . The relax- 
ation time of the free induction process corresponds to 
the lowest order (mean relaxation time) approximation, 
which provides the best single exponential approxima- 
tion, M{t) « e~*/^2 ^ by setting T^ = M(0). Thus, 



t; 



t{1 + ,) 



yi + i^TSiuy - 1] + ^[^/i + (rSiuy - 1] 



(8) 



The free induction decay of magnetization is a result 
of coherent and incoherent spin dephasing. When in- 
stead a spin echo experiment is performed, a 180° pulse 
applied at i/2 reflects the spin phase which eliminates 
the coherent contribution of spin dephasing at the echo 
time t. Time evolution of the local magnetization be- 
fore the 180° pulse is that of a free induction decay 
m{x,t/2) = exp[(DV2 + iuj{x))t/2]m{x,0). After the 
phase reflection which is described by the complex ad- 
joint, m{x,t/2) — > m*{x,t/2) there is again a free 
induction decay, i.e. the global magnetization MsE{t) 
after the echo time t is 

MsE{t) = V-^ / dV cyi])[[Dy'^ +iuj{x))t/2]x 
Jv 

eyiY>[{DV^ -iuj{x))t/2]m{x,{)) . (9) 

Again replacing DV^ by D, differentiating this relation 
and exploiting the structure of the operator D we obtain 
a general relation between the contribution of solely in- 
coherent spin dephasing on relaxation (Mse) and that 
of coherent and incoherent dephasing (M) 

dtMsEit) = ^ (I M{t/2) p -MsE{t)) , (10) 

which is solved by 

MsEit) ^'-^(r + jy^ e«/- I M(e/2) pj . (11) 

For the free induction decay being given by M(t) w 
S;y=i /i' exp(— r,yt) one obtains 



MsE{t) 



N 



^^ 1 - tR^^ 

/i,!/ — 1 ^ 



^-R,.A _ ^i? e-*/^ 



where i?^^ = (r^ + r^)/2. 

The determination of the relaxation time T2 of Mse (t) 
depends on the experimental setup. When T2 is deter- 
mined from a single echo time i or a multi-echo sequence 
(|l7| with an inter-echo time t, then T2 is usually obtained 
from 



T, 



-HMsEit))/t 



(13) 



However, since MsE{t) is not single exponential already 
in the simplest case, A'^ = 1, the above definition renders 
T2 effectively dependent on the echo time t. 

When T2 is determined from several experiments 
with varying echo times the best single exponen- 
tial approximation, MsE{t) ~ e^*'"^^, is required. 
Within the mean relaxation time approximation, T2 — 
J^dtMsEit)/MsE{0), Eq. (|lD gives 



To = T,* 



(14) 



(12) 



Equations (S) and (O) allow the analysis of various 
limiting cases. The motional narrowing regime is given 
for T^uP')^!'^ ^ 1. Then, the relaxation rate of the free in- 
duction process is I/T2* = rf^Lo^ix)) with the spatial vari- 
ance of the spin precession frequency (a;2(x)) = <;5uj'^ /2. 
This relation is the well known general result for the 
transverse relaxtion rate obtained for the motional nar- 
rowing limit. 

For very long correlation times one obtains l/T^ = 
2<,5u) and I/T2 — 0, i.e relaxation is due solely to co- 
herent spin dephasing. For very short correlation times 
one obtains Tj* « T2, i.e. in this limit the free induction 
decay is mainly due to incoherent spin dephasing. 

Another approximation is possible for small relative 
intracapillary blood volumes (<j ^ 1) and in the regime 
tSljj < 1, which is realistic in cardiac tissue. Equation 
(|) then provides l/r2* = <rr-i(^l -t- (t(5w)2 - 1), i.e. 
the rate of the free induction decay is proportional to the 
intracapillary volume. 

Kennan and Gore B determined relaxation rates from 
Monte Carlo simulations of transverse polarized spins as 
a function of the diffusion coefficient. In Figs. 1 and 2 
their results are compared to our theoretical predictions 
from Eqs. (||), ([l2|), and (p3|), using their tissue parame- 
ters and A'^ = 1, which proved to be a sufficient approx- 
imation. They demonstrate a close similarity of analyti- 
cally obtained T2 vs D and T2 vs D curves, and those of 
simulations. The correct asymptotic behaviour for very 
long and short correlation times, as derived above, is evi- 
dent in both figures. Also the location of the maximum of 
the T2 vs D curves is congruent. Figure 1 demonstrates 
furthermore that the dependence of T2 on the echo time 



is described correctly. Figure 2 exhibits the correct de- 
pendence on the capillary radius. 

We recently performed Tj mapping in hearts of vol- 
unteers p before and after application of the coronary 
dilator dipyridamol, which increases coronary flow by a 
factor of about 5 without affecting oxygen consumption. 
This decreases the deoxyhemoglobin content and one ob- 
serves an increase of T2* w 36ms by 17%. The theory 
presented above predicts an increase of 13 to 18%, which 
is in good agreement with the experimental data. 

To summarize, we developed an analytical theory for 
the description of spin dephasing in the inhomogcncous 
magnetic field of the capillary network of myocardium. 
The underlying strong collision approach is justified in 
the cardiac muscle and, most probably, albeit with mod- 
ifications according to different vascular architectures, for 
many other tissues as well. Note that this approach is not 
only valid when spin dephasing occurs on a much longer 
time scale than the local field fluctuations, which is the 
basis for the strong coUison approximation, but also in 
the static dephasing regime. 

In our applications we considered only the single- 
exponential (mean relaxation time) approximation of the 
full relaxation curve M(t), which proved sufficient for the 
parameter regime in question. However, in particular for 
short spin echo times, it can be important to describe 
the short time behavior of M{t), e. g. ■^M{t)\t=o = 0, 
in more detail. This is readily possible using the multi- 
exponential approximation for M{t), see also |1^ . 
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FIG. 1. Relaxation rates of the free induction decay of 
a gradient echo (I/T2 , GE) and spin echo experiment (1/12, 
SE) as function of the diffusion coefficient. For comparison 
with simulation data (sim.), the tissue parameters were ob- 
tained from Ref. |8| as: relative intracapillary blood volume 
? = 5%, capillary radius Re — 2.5/im, intracapillary magne- 
tization of 1.6 mgauss [Slu = 269rad s~^). The theoretical 
results (solid lines) were obtained from Eq. (H) for the free 
induction decay,and from Eqs.([L2|) and (113) for the spin echo. 
Two echo times (TE) were analyzed. 
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FIG. 2. Relaxation rates of free induction decay {I/T2, 
GE) and spin echo experiments {I/T2, echo tiine= 50 ms, 
SE). Simulation (sim.) and theoretical (solid lines)results are 
compared for a capillary radius of 2.5/im and 7.5pm. Other 
parameters are as in Fig. 1. 



